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We have studied the temperature and magnetic field dependence of the electrical resistance of 
mesoscopic, tens of nanometers thick multigraphene samples as a function of a bias voltage applied 
perpendicular to the graphene planes. We found that the resistance changes asymmetrically with 
the bias voltage sign. For large and negative bias voltages the resistance shows a non-percolative 
superconducting-like transition at T ~ 15 ... 20 K. A large enough magnetic field suppresses the 
transition. 
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I. INTRODUCTION 

Recently published studies show that electrostatic car- 
rier accumulation is an interesting tool to trigger new 
states of matter at certain interfaces. Electric field in- 
duced quasi-two dimensional (2D) superconductivity has 
been reported in the insulating oxide SrTi03ji on an 
atomically flat ZrNCl film 2 ", on KTaC>3 single crystals^, 
in LaTiOa/SrTiOa hetero-structures^, with critical tem- 
peratures ranging from 0.05 K to 15 K, after inducing 2D 
carrier densities between ~ 10 13 cm -2 and ~ 10 14 cm~ 2 , 
upon material or interface studied. The idea to trigger 
2D superconductivity by inducing excess charge at the 
surface of a dielectric sample by the application of an 
electric field was proposed already in Ref. [H for the par- 
ticular case of the surface region of L^CuCU assuming a 
weak coupling between the Cu02 layers. 

The possible existence of superconductivity in graphite 
and in graphene has been a subject of experimental^ as 
well as theoretical interest in the last 10 years. For ex- 
ample, p-type high-T c superconductivity has been pre- 
dicted to occur in inhomogeneous regions of the graphite 
structure^ or d— wave superconductivity based on res- 
onance valence bonds&. Following a BCS approach in 
two dimensions energy gap values at K of the order of 
60 K have been obtained if the density of conduction elec- 
trons per graphene plane increases to n ~ 10 14 cm -2 ,- 
in agreement with the theoretical predictions based on 
different others approaches ! 1 ^ 11 Also room temperature 
superconductivity with a d + id pairing symmetry has 
been predicted to occur in doped graphene with a carrier 
concentration n > 10 14 cm" 2 .— 

We expect that if one can trigger superconductiv- 
ity in quasi-2D regions of graphite as well as at doped 
surfaces 13 it may exist at much higher temperatures 
than in the basically 3D superconducting intercalated 
graphitic compounds^ This expected difference in crit- 
ical temperatures and within a mean field model can be 
partially understood taking into account the role of high- 
energy phonons in the 2D graphite structure itself^ Since 
most of the theoretical predictions emphasize that super- 
conductivity should be possible under the premise that 
the carrier density per graphene area is high enough to 



reach T c > 1 K^r— it is then appealing to use the electro- 
static charge doping to increase n in graphite without dis- 
turbing its quasi-2D dimensionality We note that the in- 
trinsic carrier density of defect-free graphene layers inside 
graphite is n < 10 9 cm" 2 . 15 ' 16 However, defects and/or 
hydrogen doping within regions at interfaces 1 ^ or at the 
graphite surface may show much larger carrier density, 
e.g. n ~ 10 11 cm -2 , as literature values indicate. There- 
fore, if we can increase the carrier density above some 
critical value we may expect to trigger superconductivity 
in those regions with carrier densities just below a certain 
value. We note further that the surface region of graphite 
appears then to be especially interesting according to re- 
cently published theoretical worki£ that emphasizes that 
this may show a topologically protected flat band pro- 
moting superconductivity at very high temperatures. 



II. EXPERIMENTAL DETAILS AND SAMPLES 
CHARACTERISTICS 

We have studied five multigraphene (thin graphite 
flakes) samples named SI . . . S5. Samples SI and S3 were 
graphite flakes of ~ 40 /jm length, see Fig. [TJ and with a 
non homogeneous thickness ranging between 20 nm and 
40 nm. Samples S2 and S5 showed a homogeneous thick- 
ness of 40 nm and ~ 15 /im length. Sample S4 of 90 nm 
thickness was the thickest of all studied samples in this 
work. All samples were prepared by a rubbing method 
and fixed on the 150 nm thick SiN part of a doped Si sub- 
strate following a similar procedure as shown in Ref. [l7l . 

The contacts of the sample were prepared by elec- 
tron lithography using Pd/Au thermal evaporated elec- 
trodes, see Fig. Q] (up right). The bias voltage was applied 
through a fifth electrode. The results shown in this work 
did not depend whether we apply the bias voltage at the 
electrode at the edge, as shown in Fig. [T] (bottom), or at 
some of the other electrodes. 

The resistance was measured with AC and DC input 
currents. In case of application of a bias voltage we used 
a DC current of 1 fiA with a Keithley 6221 current source 
and the voltage measurements were done with a Keith- 
ley 2182 DC nanovoltmctcr. The AC measurements were 
performed using a LR700 AC resistance bridge from Lin- 
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level off of R(T) of sample SI below 25 K or the features 
below 50 K in samples S2 and S3 and the metallic like 
behavior of S4 are due to the contributions of the surface 
of the sample (or sample-substrate) and/or due to the 
internal interfaces in the sample, as discussed in Refs. [l?] 
and [l^. As expected, the thicker S4 sample shows a 
metallic behavior below 100 K, which is related to the 
higher number of internal interfaces as reported in detail 
m Ref.[T3. These internal interfaces have a much larger 
carrier density as the graphene layers within the defect 
free regions of graphite. Therefore, we expect that in 
thicker samples the shielding of the electric field is much 
larger than in thinner ones, as we describe below. 



FIG. 1. The upper-left picture shows an optical image of the 
sample SI just after the rubbing process on the SiN substrate. 
On the upper-right picture, the same sample after the lithog- 
raphy and Pt-Au thermal evaporation. The voltage, current 
and bias voltage electrodes are indicated as they have been 
used in the whole set of measurements. In the bottom pic- 
ture, a sketch of the configuration used in this work is shown. 
MG represents the multigraphene sample, and the labels V 
and / the voltage and current electrodes. For positive bias 
voltages V g > 0, the positive electrode is connected to one of 
the contacts on the sample surface as shown in the sketch. 
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FIG. 2. (a) Temperature dependence of the normalized re- 
sistance of four samples at zero bias voltage and zero applied 
magnetic field. The amplitude of the input AC current was 
1/iA. (b) Bias voltage dependence of the resistance for sam- 
ples SI and S3 (left y— axis), S2 and S4 (right y— axis) at 2K 
and no applied magnetic field. 



ear Research. 

The temperature dependence of the resistance R(T) at 
zero bias voltage has been measured for every sample and 
it is shown in Fig. Uta). SI, S2 and S3 samples show a 
semiconducting behavior intrinsic to the graphite Bernal 
structure according to the work done in Ref. [l9l The 



III. BIAS VOLTAGE DEPENDENCE OF THE 
RESISTANCE 

A. Electrostatic screening in multigraphene 
samples 

The electrostatic screening in multigraphene samples 
is essential since without a large amplitude of the elec- 
tric field in the sample through the applied bias voltage 
we will sec no effect on the electrical resistance. To es- 
timate the penetration depth of the applied electric field 
inside our samples we use a recently published theoreti- 
cal paper—, which emphasizes that the actual screening 
depth depends on the experimental conditions, in partic- 
ular the actual doping of the sample. The intrinsic carrier 
density in each graphene layer in the graphite structure 
without defects and interfaces should be n < 10 9 cm" 2 
at the temperature of our experiments ! 15 ' 16 ' 21 The exact 
value of n for the graphene layers in each of the mea- 
sured samples is not really well known just because the 
samples are not free from defects and impurities as, e.g., 
hydrogen. Due to the relatively low temperatures of our 
measurements we take as upper limit n ~ 10 8 cm" 2 for 
the graphene layers not involved in the internal interfaces 
or at the surface. In case the carrier density n of the non- 
defective graphene layers is smaller, then the larger will 
be the penetration depth of the electric field in the sam- 
ple. Although Ref. [20| assumed a zero energy gap for 
Bernal graphite, the magnitude of the used voltages in 
our experiments is much larger than the experimentally 
obtained semiconducting gap of ~ 40 meVJ^ Therefore 
and following Ref. [2(] we estimate that at T < 30 K the 
effective penetration depth should be equal to at least 
7 graphene layers or about 2.4 nm. That would mean 
that the electric field mainly influences the near surface 
region of the samples, in case that no internal interfaces 
with much larger carrier density exist. Otherwise, as we 
will see for sample S4, these interfaces appear to largely 
screen the applied field and no effect is observed after 
application of a bias voltage. 

Our estimate is basically in agreement with other 
theoretical work, which showed that the electric field 
should be screened within a few layers from the sample 
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surfaces^ Therefore, an estimation for the electric field 
penetration depth of at least 3 to 10 graphene layers from 
the surface appears realistic. Note that due to the large 
electrical resistance of the SiN layer and taking into ac- 
count its capacitance (~ 10 nF), of the same order as 
the effective capacitance of the multigraphene samples, 
we expect large equilibrium relaxation times in the range 
of t = RC > 5 x 10 3 s (R, C are the total resistance and 
capacitance parallel to the applied electric field). 

Finally, we note that the bias voltage is applied be- 
tween the top of the Si substrate and the top surface of 
the samples. This means that we have a distance between 
electrodes of 150 nm plus the thickness of the samples, 
i.e. between 170 nm to 240 nm. This effective thickness 
difference of < 40% between the samples, and within a 
rough linear capacitor model, would indicate a reduction 
of the same order in the electric field amplitude for the 
same bias voltage between the thinnest and thickest sam- 
ples. We will see below, however, that it is not only the 
thickness the parameter that determines the amplitude of 
the effect of the electric field but the highly conducting 
interfaces, which have a mean distance between 30 nm 
and 100 nm inside the graphite structure of the studied 
samples 



B. Bias voltage dependence at constant 
temperature 

Figure [SJb) shows the resistance vs. applied bias volt- 
age V g at a constant temperature of 2 K without apply- 
ing any magnetic field. Starting from zero bias voltage 
we increased to positive V g values, then we returned to 
zero and measured the resistance for the negative range 
of V g . Each point in Figure |5Jb) was measured after 
reaching equilibrium at each point till the resistance did 
no change with time. The curves shown in Fig. [UJb) are 
reversible, demonstrating the electrostatic nature of the 
observed effects. 

The samples SI, S2 and S3 show an asymmetric behav- 
ior respect to zero voltage. For positive voltage, i.e. the 
electric field direction points towards the sample and sub- 
strate, the resistance does not change significantly. In the 
case of sample SI it slightly decreases after +50V. The 
resistance of S2 does not change up to +60V, and finally 
the S3 resistance increases slightly for V g > +50 V. 

For negative V g , however, a remarkable decrease of the 
resistance appears. A minimum value in the measured 
resistance is reached for slightly different negative bias 
voltages upon sample. For sample S2 the minimum of 
resistance appears at -50V, for S3 at -60V and for SI, 
which shows the largest change, at ~ —75 V. As will 
become clear in the discussions below, these differences, 
as well as that in the relative decrease of resistance at 
a given bias voltage, are related to the carrier inhomo- 
gencity of the multigraphene samples. Note that for the 
thickest sample S4 the resistance does not show any sig- 
nificative change with V g . This is a clear indication for 
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FIG. 3. Temperature dependence of the resistance for four of 
the multigraphene samples at V g = — 60 V and at no applied 
magnetic field. 



the screening effect of the electric field by the internal 
interfaces mentioned previously. Would the electric field 
amplitude decreases simply proportional to the sample 
thickness, we expect to see a decrease in the resistance 
within the scanned V g range. A similar behavior of the 
resistance of the multigraphene samples SI to S4 with 
the applied bias voltage was partially reported in Ref. l23l . 
However, in that work no results on the temperature or 
the magnetic field dependence of the resistance under a 
bias voltage were reported. 



C. Temperature dependence of the resistance 
under a bias voltage 

The results of the resistance vs. bias voltage shown in 
Fig. [2] were obtained at a constant temperature of 2K. 
Taken into account these results, the aim of the next ex- 
periment was to investigate the change of the resistance 
with temperature at constant applied bias voltage. Fig- 
ure[3]shows the temperature dependence of the resistance 
of four samples without and under V g = —60 V. As one 
can easily realize and as it was shown before in Fig. [2th), 
the resistance is in general smaller under -60 V. But re- 
markable is the fact that in all the samples a clear step 
like transition below 20 K appears. The details of the 
transition depends on the selected sample, again a sign 
of the existence of inhomogencities in the samples. 

Furthermore, the effect of the gate voltage values on 
the temperature dependence of the resistance was inves- 
tigated in detail. Figures H]and [5] show the results mea- 
sured for samples SI and S2, respectively. In the case of 
sample SI, the temperature behavior of the value as well 
as the temperature dependence of the resistance R(T) 
does not change significantly for 
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-40 V a small dip in R(T) appears at T ~ 17 K. 
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FIG. 4. Temperature dependence of the resistance of sample 
SI at different constant values of V g and no applied magnetic 
field. 



For V g < —40V the small dip develops in a clear step 
at similar temperatures with less than 2 K transition 
width. The overall shape of these curves suggests the ex- 
istence of a non-percolative superconducting transition, 
as the measured behavior with magnetic field also indi- 
cates. The curve obtained at V g = —60 V with a drop 
at 15 K and an upturn below it are clear indications of a 
non-uniform channel with superconducting regions con- 
nected in series with normal ones in which partial charge 
localization effects occur. The clear asymmetric behavior 
of the resistance with V g is recognized also in Fig. fJJ 

The lower the bias voltage the clearer the transition 
triggered by the change in the carrier density. Note that 
the apparent transition temperature does not shift sig- 
nificantly with the applied bias voltage. The overall be- 
havior of the resistance R(T) with applied V g strongly 
resembles in several details the one obtained for ZrNCl 
flat thin film in Ref. [U, even in the obtained transition 
temperature. However, with one important difference, 
namely that in the multigraphene samples no supercon- 
ducting percolating path is reached. We ascribe this be- 
havior to the carrier inhomogeneities of the samples^, a 
fact related to the large sensitivity of the carrier density 
of the graphite structure to defects and impurities 

A similar, but not identical, behavior is obtained for 
sample S2, see Fig. [5] This sample shows already at 
—20 V applied bias voltage a clear step in the resistance 
starting at ~ 17 K. At higher negative voltages, how- 
ever, the resistance starts to show signs of a transition 
at higher temperatures. It is clear that the details of the 
effects observed by the applied electric field depend on 
the multigraphene sample, again an evidence for the role 
that inhomogeneities play in the observed phenomenon. 
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FIG. 5. Temperature dependence of the resistance of sample 
S2 at different constant values of V g and no applied magnetic 
field. 



D. Temperature dependence of the resistance 
under a magnetic field and bias voltage 

To characterize the magnetic field behavior of the ob- 
served transition as a function of temperature we fixed 
the bias voltage to a large enough value taking into ac- 
count the results shown in the last sections. In what 
follows we discuss mainly the results of sample SI for 
both, magnetic field applied normal and parallel to the 
graphene planes. Results for the other samples are qual- 
itatively similar. Figure [5] shows the dependence of the 
resistance (absolute (a) and relative (b)) with tempera- 
ture at different applied fields normal to the main area 
of the sample and at a fixed bias voltage of -100 V. The 
measurements were performed as follows. At 2K a bias 
voltage of —100 V was applied. After reaching equilib- 
rium the corresponding magnetic field was applied as well 
and the resistance was measured to 25 K. At 25 K the 
magnetic field was removed, the sample was cooled again 
to 2 K again and we repeated the same procedure for 
others magnetic field values. The bias voltage has been 
applied all the time. 

As expected for this field direction, the background re- 
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FIG. 6. Temperature dependence of the resistance of sample FIG. 7. Temperature dependence of the resistance of sample 
SI at different constant values of magnetic field applied nor- SI at different constant values of the field applied parallel to 
mal to the graphene planes and at V g = -100 V. (a) Absolute the graphene planes and at V g = -100 V. 
values of the resistance, (b) Normalized values of the resis- 
tance by the resistance at 25 K vs. temperature at the same 
applied fields as in (a). 



sistancc increases with field due to the usual magnctorc- 
sistance of the multigraphene samples, see Fig. [6{a). As 
shown in Fig. [6] a field of 0.2 T is enough to suppress com- 
pletely the transition at 15 K. This suppression remains 
to a field of 1.5 T. At a field of 3 T and higher the transi- 
tion appears again at the same temperature but slightly 
broader and nearly vanishing at 8 T, see Fig. |6jb). This 
non monotonous field behavior is observed for fields nor- 
mal to the graphene planes only. Note that the temper- 
ature of the transition does not change significantly with 
applied field. The transition is rather unconventional be- 
cause the magnetic field affects mainly the relative step 
height of the transition, see Fig. EKb). 

Figure [JJ shows the resistance vs. temperature at 
V g = —100 V and at different magnetic fields applied par- 
allel to the graphene planes of the sample. The misalign- 
ment of the field is less than 0.5°. The measured data 
show that the transition remains unaffected by a paral- 
lel field of 0.2 T, in contrast to the normal field result. 
At higher fields, however, it is suppressed monotonously 
without any sign for a reentrance like in the other field 
direction. In Fig. [8] we compare the results for both field 
direction by plotting the minimum resistance just below 
the transition normalized by the resistance value at 25 K 
vs. applied field. As for the normal applied fields, for 
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parallel applied fields the transition does not shift signif- 
icantly to lower temperatures and the resistance shows a 
minimum just below the transition. Note that the resis- 
tance above ~ 17 K does not change practically with field 
in agreement with the fact that the magnetoresistance of 
graphite depends mainly on the normal field component 
to the graphenc planes^ 



IV. DISCUSSION 

In this section we discuss the rather unconventional 
behavior of the electrical resistance under a bias voltage 
and its possible origin. We divide the discussion in three 
parts. In the first part (A) we show experimental ev- 
idence that supports an increase in the carrier density 
after application of a bias voltage in our multigraphene 
samples. In part (B) we discuss to what extent the ob- 
served transition can be considered as due to supercon- 
ductivity or other kind of electronic transition. We em- 
phasize the resemblance of the observed behavior with 
that obtained in the past for granular superconducting 
thin metallic films. We shortly review the main results of 
recently published work on graphite samples that suggest 
the existence of granular superconductivity in this mate- 
rial. In part (C) of this section we speculate on possible 
explanations for the anomalous magnetic field behavior 
of the transition, assuming that it is indeed related to 
superconductivity. 



A. Increase of carrier density with bias voltage 

The transition observed at T < 20 K in different 
multigraphene samples after the application of a negative 
bias voltage indicates that this happens after increasing 
enough the carrier density and the Fermi level. Taking 
into account the expected field penetration depth, we 
expect that the superconducting path is located at the 
upper near surface region of the samples. A way to check 
whether there is a real increase in the carrier density with 
applied bias voltage is to measure the Shubnikov-dc Haas 
(SdH) oscillations of the magnetoresistance. Figure [9ja) 
shows the magnetoresistance of sample SI at 2 K with- 
out and with an applied voltage. From these data we 
plot in Fig. HJb) the first field derivative of the magne- 
toresistance vs. inverse field. At no applied bias voltage 
the data reveal no SdH oscillations in the shown field 
range and within experimental error, in agreement with 
the rather small carrier density of the graphene planes 
n < 10 9 cnrlMS 

At an applied bias voltage of -100 V clear SdH oscil- 
lations are observed in the first field derivative for fields 
above ~ 4 T. From the obtained period of the oscilla- 
tions (lowest frequency / ~ 46 T) we estimate a 2D car- 
rier density n ~ 2.2 x 10 12 cm -2 , an order of magnitude 
larger that the one obtained for the bulk graphite sam- 
ple (continuous line in Fig. E^b)) from which the multi- 
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FIG. 9. (a) Magnetoresistance of sample SI for normal to the 
graphene planes field, at 2 K and without and with an applied 
bias voltage, (b) First derivative of the magnetoresistance vs. 
inverse field from the data of (a) and for a bulk graphite 
sample of size 2 x 1 x 0.2 mm 3 at 4 K (upper-right axes). 



graphene samples were prepared. We stress that the car- 
rier density obtained for the bulk sample is not intrinsic 
of the Bcrnal graphite structure but it originates mainly 
at the internal interfaces (or other defective regions, see 
Ref. [Tol ) commonly found in the used highly oriented py- 
rolytic graphite samples.— As shown recently^ the SdH 
oscillations in multigraphene samples depend on the sam- 
ple position indicating the existence of inhomogeneities in 
the carrier concentration within micrometers. This pic- 
ture agrees with the results obtained from electric field 
microscopy (EFM) that revealed submicrometer domain- 
like carrier density distributions in graphite surfaces ^ 
On going experiments using an array of electrodes dis- 
tributed on the same sample provide further evidence for 
the inhomogeneous response of the multigraphene sam- 
ples after application of the bias voltage. This result also 
indicates the existence of inhomogeneities in the carrier 
distribution within the sample. 



B. Can the observed transition be related to 
superconductivity? 

The obvious question that needs an answer is the origin 
of the observed transition in the resistance when the car- 
rier density increases. Leaving apart superconductivity, 
it appears difficult to find the origin in a kind of electron 
derealization decreasing temperature. Also the fact that 
the magnetic field removes the transition appears to be 
at odd with electron localization effects. Note that a par- 
allel magnetic field does not change the resistance above 
the transition, see Fig. [Jj indicating clearly the absence of 
any Lorenz-force driven effect or a change in the electron 
system at all above the critical temperature. 

One may speculate that the observed transition might 
be related to an electric field induced sliding motion of a 
charge-density wave (CDW), which may lower the resis- 
tance, and the existence of a strong coupling between a 
magnetic field and the CDW structure, as has been ob- 
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served in NbSe3, for exampl o 26 ' . However, in the case 
we present here two observations appear to be at odd 
with a CDW origin. Namely, the transition at T ~ 15 K 
is sharp and within less than 3 K in contrast to the very 
broad decrease of the resistance due to the move of a 
CDW decreasing temperature under an applied electric 
field 2 ^. Second, we did not find any signature of non- 
linearity of the resistance with current I below 10 fiA 
at any applied magnetic or electric field (currents above 
this value can destroy the sample and were not applied) , 
in contrast to the observations in CDW transition o 27 ' 28 . 
Note that in case of a CDW the resistance becomes non- 
linear at progressively lower currents as the magnetic field 
is increased 2 ^. 

The absence of nonlinearity at low currents may appear 
to be also at odd with the existence of granular supercon- 
ductivity. However, the appearance of nonlinearities in 
the resistance in this case depends actually on the sample 
granular characteristics and current range. For compar- 
ison we refer to the extended work done on the onset of 
superconductivity in ultrathin granular metal films, see 
for example Refs. and HO- The similarity of our results 
of the resistance vs. temperature for constant applied 
electric fields (see Figs. 3 and 4) with those found in 
non-percolative granular superconducting films^ 9 - is ap- 
pealing. We note that in the case of ultrathin super- 
conducting Ga films the I — V characteristics are rather 
linear at low enough currents and in films where no global 
superconductivity occurs 2 ^. 

On the other hand the existence of granular high- 
temperature superconductivity in regions within inter- 
nal interfaces found in highly oriented pyrolytic graphite, 
has been proposed recently^ 7 - as a possible answer to a 
series of experimental hints on its existence in oriented 
graphite samples published the last years. In short, the 
main experimental observations that suggest the exis- 
tence of non-percolative superconducting regions in ori- 
ented graphite samples are: (1) Magnetic field driven 
metal-insulator transition (MIT) 31 i 32 . We should note 
that the MIT in graphite has been also interpreted in 
terms of the two band model using several magnetic field 
and temperature dependent free parameters 3 - 3 -. How- 
ever, this "transition" does not exist when the oriented 
graphite sample gets thinner than ~ 70 nm and therefore 
its origin is not intrinsic of the graphite ideal structure 1 ^; 
(2) Superconducting-like hysteresis loops in the mag- 
netization of pure HOPG samples^; (3) Anomalous 
hysteresis in the low-field magnetoresistance of meso- 
scopic graphite samples, similar to those found in granu- 
lar superconductors 3 ^^; (4) Superconducting-like mag- 
netic field driven transition in the resistance when mea- 
sured near or at an internal interface of a graphite flake 37 ; 
(5) Anomalous temperature dependence of the zero field 
resistance of oriented graphite samples, which follows 
an exponential, thermally activated behavio r 19 ! 38 simi- 
lar to the one found in granular superconductors as in 
granular Al-Ge 3 ^, for example. (6) Finally, Josephson- 
like I — V characteristic curves found in especially pre- 



pared thin graphite transmission electron microscope 
(TEM) lamellas 4 ^. Apart from the hints in pure ori- 
ented graphite samples, experimental evidence for high- 
temperature superconductivity has been also reported in 
composite graphite - sulphur powders4i~— All this evi- 
dence indicates clearly an intrinsic tendency of the phe- 
nomenon of superconductivity to appear in certain re- 
gions of graphite after doping. 

In order to reach a state of zero resistance in some part 
of the graphite sample, it is clear that one needs to have 
a doped region large enough that one or several Joseph- 
son coupled superconducting patches exist between the 
voltage electrodes. This appears still difficult to reach 
applying only an electric field in the samples used in this 
work. On the other hand results on thin TEM lamel- 
lae prepared with a large number of internal interfaces, 40 
which show a carrier density similar or higher than the 
one we obtained in our samples after applying an elec- 
tric field, shows a current dependent zero resistance state 
below 20 K. Taking into account this evidence and from 
the comparison with literature it appears natural to in- 
terpret the observed transition with electric field in our 
multigraphene samples as a sign of granular supercon- 
ductivity. 



C. Anomalous magnetic field behavior of the 
transition 

The magnetic field behavior of the transition for both 
field directions, i.e. the resistance below the transition 
increases with field without a clear decrease of the tran- 
sition temperature within the used field range, suggests 
that a filamentary superconducting path produced by the 
applied bias voltage, is affected by the magnetic field. 
Future experiments should clarify whether the super- 
conducting channel connects localized superconducting 
"grains" or patches with much higher critical tempera- 
tures and fields. 

The differences in the behavior observed as a function 
of field direction are also peculiar. From one side or- 
bital effects appear to be important taking into account 
the suppression of the transition for a field of 0.2 T (for 
sample SI, for example) applied normal to the graphene 
planes only. Also the reentrance of the transition ob- 
served only for fields normal to the graphene planes ap- 
pears to have an orbital character. We note that a similar 
effect has been observed through the measurement of the 
conductance of a high-mobility 2D electron gas between 
superconducting contacts at high fields applied normal 
to the main 2D area4^ This reentrance or increase in 
the conductance with magnetic field was explained argu- 
ing the increase in the probability of Andreev reflections 
above a certain field4^ A comparison of our results with 
those from Ref . is permissible because the carriers mo- 
bility in the graphene layers of our samples is huge 1 ^ and 
also hints for the influence of Andreev reflections in the 
magnetoresistance have been recently reported^ 
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On the other hand we should take into account that 
the transition is also affected by a parallel field. The 
absolute difference at low fields (for sample SI) between 
the two field directions cannot be explained by the small 
misalignment in the applied field direction and indicates 
that other mechanisms than orbital may play a role in the 
triggered superconducting state. The difference in the 
resistance behavior as a function of the field direction 
might be understood assuming the interaction between 
preexisting ferromagnetic moment and spin-polarized su- 
perconducting currents (SC) as was conjectured for su- 
perconducting graphite doped with sulfur4£ The spin- 
polarized SC currents, associated with either p-wave 
spin-triplet 7 - or <i-wave superconducting order parametei— 
with breaking time-reversal symmetry, exert a torque on 
the ferromagnetic moments leading to their rotatio n 47 i 48 
towards the graphene planes. Then, strong enough in- 
plane fields (B > 0.2 T) orients the ferromagnetic mo- 
ments and suppresses the Josephson tunneling between 
superconducting "grains" , leading to the minimum in the 
resistance and the reentrance to the normal state at lower 
temperatures, see Fig. [7] The complete suppression of 
the superconducting transition by the out-of-plane mag- 
netic field of 0.2 T as well as the reappearance of the 
superconducting transition in (quantized) magnetic fields 
should be taken as a clear manifestation of orbital effects, 
absent in the parallel field geometry. Noting also that 
high enough perpendicular magnetic fields can rotate the 
ferromagnetic moments in the field direction suppressing 
the Josephson "intergrain" coupling, both orbital and 
"spin" effects play a role for this field direction. Evi- 
dence for a ferromagnetic response can be obtained from 
the magnetoresistance for parallel field direction, as was 
published recently for a similar multigraphene sampled 



V. CONCLUSION 

We studied the behavior of the resistance of several 
multigraphene samples as a function of temperature and 
magnetic field and under the influence of a bias voltage 
applied normal to the graphene planes. Taking into 
account relevant literature on granular superconductors 
as well as the one obtained recently for the internal in- 
terfaces in graphite samples, it appears natural to as- 
sume that the transition in the resistance that develops 
at T ~ 15 K with negative bias voltage is related to a 
non-percolative superconducting-like state. The results 
indicate that a negative voltage increases the carrier den- 
sity and the Fermi level in certain regions of the samples 
achieving the necessary minimum density to trigger su- 
perconductivity. Shubnikov-de Haas oscillations in the 
magnetoresistance and above 4 T indicate an increase 
in average of the carrier density of one order of magni- 
tude relative to the one obtained in bulk graphite sam- 
ples (without bias voltage). A large enough magnetic 
field in both directions suppresses the transition. How- 
ever, we found quantitative and qualitative differences in 
the behavior of the resistance upon field direction. The 
overall results also indicate that the carrier density is 
not homogeneous in the studied multigraphene samples 
in agreement with previous reports based on EFM and 
SdH oscillations experiments. 
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